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SYMMETRY  IN  ELECTROMAGNETIC  SCATTERING 
FROM  A  CLUSTER  OF  SPHERES 


L  INTRODUCTION. 

This  paper  is  a  sequel  to  one  we  published  in  Applied  Optics.1  Henceforth,  the 
previous  article  will  be  referred  to  as  Part  I.  We  shall  briefly  describe  the  subject  as  it  was  in 
Part  I  and  then  extend  it  here. 

As  stated  in  Part  I,  we  are  investigating  the  problem  of  electromagnetic  scattering  by 
particles  of  arbitrary  shape.  We  have  built  a  model  scatterer,  assuming  for  computational  simplicity, 
to  know  the  complex  refractive  indices  and  radii  of  the  cluster  of  homogeneous  nonmagnetic 
spheres  of  which  it  consists.  The  geometry  of  the  cluster  is  rigid  and  we  assume  that  the  incident 
wave  is  a  circularly  polarized  plane  wave  of  wave-vector  k  and  frequency  oj.  Such  an  arrange¬ 
ment,  where  in  the  cluster  the  spheres  are  homogeneous  and  nonmagnetic,  allowed  us  to  transform 
the  vector  problem  (the  scattering  problem  has  an  intrinsic  vector  nature)  into  a  pair  of  scalar 
ones  using  the  Debye  potentials  to  describe  the  electromagnetic  field.  In  this  way,  we  have  simpli¬ 
fied  the  problem  and  have  been  able  to  compute  the  absorption,  scattering,  and  total  cross  sections. 

The  computations  of  the  cross  sections  of  the  cluster  and,  therefore,  the  convergence 
of  the  scattered  fields  require  solution  of  rather  big  systems  whose  order  is  given  by  N(£max  +  l)2, 
N  being  the  number  of  the  spheres  in  the  cluster,  and  £„,„  the  highest  value  of  £  that  is  included 
in  the  expansions  of  Debye  potentials.  If,  however,  the  cluster  has  a  symmetry  group,  the  above 
systems  can  be  factorized  through  group-theory  techniques  similar  to  those  used  in  quantum- 
mechanical  problems.  It  is  shown  that  the  application  of  group  theory  to  scattering  produces  a 
nonhomogeneous  system  of  equations  rather  than  secular  determinants.  These  nonhomogeneous 
equations  more  closely  represent  real-world  scattering  by  particles  than  is  generally  obtained. 

II.  THE  SCATTERING  EQUATIONS  FOR  GENERAL  DIRECTION  OF  INCIDENCE. 

In  Part  I,  it  was  assumed  that  the  incident  wave  propagated  along  the  positiv  z-axis, 
the  direction  of  incidence  being  then  changed  by  rotating  the  cluster  as  a  whole.  For  our  present 
purposes,  we  shall  refer  the  cluster  to  a  fixed  system  of  axes  and  choose  the  direction  of  incidence 
through  the  direction  cosines  of  k,  the  propagation  vector  of  the  wave.  A  straightforward  calcula¬ 
tion  along  the  lines  sketched  by  Jackson2  shows  that  for  a  circularly  polarized  plane  wave  of 
wave-vector  k: 


£*(1)  “  (fi  ±  ifi>)  exp  (ik  •  I)  * 


4x£ifi(fii  ±  ifia)  '  Xftm<£)  I  *8*(l)i£(kr)  ±  r  Ax  *tm(i)it<kr) 

£m  L 


(la) 


I*(I)  -  «£*(r) 


(lb) 
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where  £(  and  &  are  unit  vectors  orthogonal  to  fc  and  to  each  other  and  (£)  and  Xgn,(i)  are 
vector  spherical  harmonics  defined  as: 


<*> 


-it *  V 

Vw  + 1) 


•  ■iwt) 


-ik  X  y 

£(B  +  1) 


(2) 


The  Debye  potentials  for  the  fields  (1)  are  then  (Newton3) 

♦St!)  -  f  g  i*d,Q5„<£>j,<kr)Yta<£)  (3.) 

<(!)-*  *J(X)  OW 

with 

dj  *  4*  [«(«  +  1)J**  ,  Qtnd)  -  (fi.  ±  fc»)  *  *fa,(k) 

When  the  above  expressions  for  p*  and  p*  are  introduced  into  the  formalism  developed  in  Part  I, 
the  final  systems  of  equations  which  determine  the  coefficients  Bj<m>  become 

[4}'' ♦  <,•„■}  <.•  ■  (4*> 

EEjVn-Wft.]  1  +  H,®  <4» 

P 

where  s£,  Tj  and  are  still  defined  as  before. 

III.  SYMMETR1ZATION  OF  DEBYE  POTENTIALS. 

Hie  symmetrization  of  the  systems  (4)  requires  us  to  symmetrize  only  Pn  and  pa, 
the  Debye  potentials  in  the  intersphere  region.  The  expansion  coefficients  of  P|  and  pj,  the 
potentials  within  the  spheres,  cancel  out  of  both  sides  of  the  resulting  equations.  It  is  convenient 
to  partition  the  cluster  in  subsets  of  spheres  which  are  related  to  each  other  by  the  symmetry 
operations  and  to  indicate  the  vector  position  of  the  ot-th  sphere  in  the  o-th  subset  by  R£,  while 
r®  *  r  -  R®.  The  expansion  of  in  terms  of  symmetry-adapted  combinations  of  Helmholtz 
solid  harmonics  is: 

*»“>***  ££  <6> 
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where  we  use  the  definitions 


<rw£«,*p>ML(:; 


(7) 


for  the  combinations  of  solid  harmonics  centered  at  the  various  sites  of  the  o-th  subset  and 


<?<£>  -V  bM«i'<k,)''tm  ® 


(8) 


for  the  combinations  of  functions  centered  at  the  origin. 


In  the  preceding  equations,  the  superscripts  p,  p  refer  to  the  p-th  row  of  the  p-th  irreducible  repre¬ 
sentation.  Note  that  ■  H^|p,  J^JP  *  J^|p  so  that  a  given  irreducible  representation  may 
appear  more  than  once  for  a  given  S.  The  symmetrization  coefficients,  ajjf™  bjj|m,  are 
easily  obtained  by  standard  techniques4, 5  through  the  use  of  the  well-known  transformation  for¬ 
mulas  for  spherical  harmonics.6  A  quite  analogous  expansion  can  be  written  for  the  coef- 
fients  being  indicated  by  °.  A  comparison  of  the  symmetrized  expansions  with  their  unsym¬ 
metrized  counterparts  yields  the  equations 

uom  ok  turn  xm  uotn  nx  nxm 

*  w  r  p 


m 


(9) 


which  allow  us  to  calculate  the  cross  sections  whose  expressions  involve  the  A®a ’s  and  the  ’s. 

xm  Xm 


IV.  FACTORIZATION  OF  THE  SYSTEMS  FOR  THE  COEFFICIENTS. 


The 

0  the  well-known 


equation  for  the  coefficients  Ajff  0  and  are  obtained  by  imposi 

n  boundary  conditions  (Newton*)  on  the  surface  of  each  sphere.  To 


imposing  on  0  and 
To  do  this,  we 


first  put  0q  and  in  a  form  involving  only  Helmholtz  solid  harmonics  centered  at  a  single  site, 
say  r£,  through  the  same  addition  formulas7, 8  we  have  already  used  for  the  unsymmetrized  case. 


The  result  for  0jj  is 

£E 

where  we  put 


ppo  ppca.(l) 
n'£  Vsm  £ 


W 


(10) 


mi 
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and: 


C.«w  (-<)'  Yf.  »»(-^(«:) 


When  we  impose  the  required  boundary  conditions  on  d  at  the  surface  of  the  ar-th  sphere  in  the 
o-th  subset  we  get  for  each  n,  p  an  equation  of  the  form: 


.  .  .  .  Tot  l-1  .  _oo.  fir  I  uprfi  upr  = 

or  afi  ££#  mm  l  *J  $m;  J  n*£W  n'£f 


if  ofW 

We  notice  that  Sj  should  be  written  as  SgTbut,  since  it  is  actually  independent  of  the  site  within 
each  subset,  the  superscript  0  has  accordingly  been  dropped.  Now,  if  equation  1 2  is  multiplied  by 
and  summed  over  a  and  m  and  fixing  the  superscripts  p,  p,  and  dropping  them  from  the 
equanon,  we  get 


nY  r 


6 

or  £  £ 


5nn'  [S£'l  +  C£n;  £'n'  |  ^n'£'  *  ~2£V  ^n; ' 


Equation  13  is  obtained  by  applying  the  following  techniques: 


ror 

£n;£V 


afi  mm 


(oa  \*  raa,  fir  rfi 
l  aim)  £m;£'m'  Vfi 


£*ir  “  (an£m)  C£m;  8V  bfi8W 

a  am 

and  where  we  have  taken  into  account  that,  if  the  irreducible  representations  are  chosen  properly 
(Altmann  and  Bradley9  and  Diamond6),  the  symmetrization  coefficients  satisfy 


E(*ry< 


,"Paa  »  6 

n'£m  nn' 


so  that 


***  *«*'  (*»*«)  [^]  *  *nYm'  "  Kr  *££'  [S£|  * 


Equations  quite  analogous  to  equation  13  hold  true  for  the 

We  recall  that  (although  never  explicitly  indicated)  the  coefficients  and  ° 
depend  on  both  the  polarization  and  on  the  direction  cosines  of  L  Furthermore,  once  then’s  and 
B' s  have  been  computed,  the  C"s  and  D's,  the  expansion  coefficients  of  the  symmetrized  Debye 
potentials  within  the  spheres,  are  calculated  through  equation  16  of  Part  I  which  are  still  valid  on 
account  of  equations  14  and  9. 


The  right-hand  side  of  equation  13  involves  the  quantities  (?gn.  which  are  easily 
identified  as  the  symmetrized  matrix  elements,  in  the  site  and  angular  momentum  representation, 
of  the  free-space  propagator  for  plane  waves.1®  Doing  this  not  only  simplifies  the  expression  of 
the  incident  wave  in  symmetrized  form  but  also  allows  the  calculation  of  the  cross  section  without 
any  approximation  other  than  the  truncation  of  the  1 -expansion  of  Debye  potentials.  Indeed  the 
scattered  potentials  are  given  by: 


2k 


aa  £m 


?  a  „r  y  Cry.1’  ftr”) 


EE 


.£  .  .aa 
1  £  fim  Y« 


a  a  £m 


(15b) 


and,  as  the  potentials  of  the  incident  wave  involve  only  Helmholtz  harmonics  centered  at  the 
origin,  it  is  convenient  to  use  again  the  addition  formulas  of  Nozawa7  to  transform  equation  1 S 
into  the  form: 


*.V  ’  -  a 


j £j  »«« 

1  d£®£*n 


aa  £m  £W 


*.<*>  =  *  2k 


aa  ftm  £#m 


(^Y*,  ^hi,)<kr>  (l6a> 

; i,d,c8c«w(<)YtaA»jV)  «« 


We  notice,  incidentally,  that  when  kRft«  1,  equation  16  reduces  to  equation  17  of  Part  I  for 
on  account  of  the  small-argument  behavior  of  the  spherical  Bessel  functions 


lim  Sa°  ,  , 

kR°«o  ®n,’*,n 

Q 


*££'  6mm' 


By  defining 
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*ao  aa 


11 


:  V :  *fc.-  «• ;  cn»h  ^  *y*.  t- .  • . 


we  get  for  the  cross  sections: 


vi.  DISCUSSION- 

The  procedure  described  in  the  preceding  sections  does  not  differ  from  that  customarily 
used  to  factorize  secular  determinants,  except  that  the  inhomogeneity  of  our  equations  forces  us  to 
solve  all  the  systems  arising  from  equation  4.  Indeed,  the  quantities  Sg>,  T^.,  on  the  left-hand 
jjide  of  equation  4  are  the  matrix  elements,  in  the  site  and  angular  momentum  representation,  of  the 
t  +  G,  the  ^-operator  accounting  for  the  scattering  power  of  a  single  sphere11'  and  G  being  the  free 
space  propagator  for  outgoing  spherical  waves.12  As  the  operator  ^  +  G  is  invariant  under  the 
symmetry  group  of  the  cluster,  its  symmetrized  matrix  elements,  i.e.,  the  quantities  within  braces 
in  equation  13,  are  actually  independent  of  the  row  index,  p.  On  the  other  hand,  the  potentials 
of  the  incident  wave  on  the  right-hand  side  of  equation  4  are  not  invariant  under  the  symmetry 
group.  According  to  general  theorem,  they  have  been  decomposed  into  parts  belonging  to  the 
rows  of  the  irreducible  representations,  but  these  symmetrized  parts  do  actually  depend  on  the 
row  index.  However,  because  of  the  above-discussed  independence  of  p  on  the  left-hand  side  of 
equation  13,  very  little  extra  computational  work  is  required  to  solve  all  systems  belonging  to  a 
multidimensional  representation.  From  a  computational  point  of  view,  most  of  the  computing 
time  is,  in  fact,  spent  to  invert  the  matrix  of  ^  +  G  and,  as  the  required  time  grows  more  quickly 
than  the  order  of  the  matrix,  the  factorization  procedure  is  always  very  useful. 

As  an  example,  let  us  consider  a  tetrahedral  cluster  with  CH4  structure  (symmetry 
group  Tj).  By  including  terms  up  to  fimax  *  2,  we  get  two  unsymmetrized  systems  of  order  45. 
However,  by  using  symmetrized  expansions,  each  one  of  the  above  systems  splits  into  one  A2 
system  of  order  3,  one  Aj  of  order  1,  two  E  of  order  4,  three  Fj  of  order  6,  and  three  F2  of 
order  5.  The  usefulness  of  the  factorization  through  group  theory  needs  no  further  comment. 

VII.  SUMMARY. 

The  work  reported  here  and  in  our  previous  papers  has  led  to  a  unique  approach  to 
calculating  the  optical  properties  of  clusters  of  molecules.  It  now  appears  possible  to  compute 
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extinction  of  a  particle  knowing  only  its  chemical  composition.  This  is  a  radical  departure  from 
the  Mie  theory,  yet  the  solution  converges  to  Mie  results  for  a  single  sphere.  The  theory  developed 
here  will  now  be  tested  using  the  X  a  -  <o,  scattered  wave  computer  program  developed  by  Johnson, 
et  at.  at  IBM  and  modified  by  the  authors  for  this  study.  The  initial  material  to  be  investigated 
will  be  sulfuric  acid,  a  material  which  is  optically  well  characterized.  This  validation  effort  will 
be  the  subject  of  a  subsequent  technical  report. 
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Box  46,  APO  New  York  09710 

Commander 

US  Army  Science  A  Technology  Center-Far  East  Office 
ATTN:  MAI  Borges  1 

APO  San  Francisco  96328 

Commander 
2d  Infantry  Division 

ATTN:  EAIDCOM  1 

APO  San  Francisco  96224 

Commander 

5th  InUntry  Division  (Mech) 

ATTN:  Division  Chemical  Officer  1 

Fort  Polk,  LA  71459 


OFFICE  OF  THE  SURGEON  GENERAL 
Headquarters 

US  Army  Medical  Research  A  Development  Command 
ATTN:  SGRD-PL  l 

Fort  Detrick,  MD  21701 

Commander 

USA  Biomedical  Laboratory 

ATTN:  SGRD-UV-L  1 

Aberdeen  Proving  Ground,  MD  21010 

US  ARMY  HEALTH  SERVICE  COMMAND 

Superintendent 
Academy  of  Health  Sciences 
US  Army 

ATTN:  HSA-CDH  1 

ATTN:  HSA-IPM  I 

Fort  Sam  Houston,  TX  78234 

US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Project  Manager  Smoke/Obecurants 

ATTN:  DRCPM-SMK  t 

Aberdeen  Proving  Ground,  MD  21005 
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DISTRIBUTION  LIST  S  (Could) 


Names 

Copies 

Names 

Copies 

Commander 

Dfaector 

US  Army  Materiel  Development  and  TtssiUnsei  Coma 

and 

RaUstic  Research  Laboratory.  ARRADCOM 

ATTN:  DRCLDC 

1 

ATTN:  DRDAR*TSB-S 

1 

ATTN:  DRCSF-P 

1 

Abardeen  Proving  Ground,  MD  21005 

5001  Eisenhower  Ave 

Alexandria,  VA  22333 

US  ARMY  ARMAMENT  MATERIEL  READINESS  COMMAND 

Commander 

Commander 

US  Army  Foreign  Science  4  Technology  Center 

US  Army  Armament  Materiel  Reedfawss  Commend 

ATTN:  DRXST-MT3 

1 

ATTN:  DRSAR-ASN 

1 

220  Seventh  St,  NE 

ATTN:  DRSAR-PDM 

1 

Charlottesville,  VA  22901 

ATTN:  DRSAR-SF 

1 

Rock  Island,  1L  61299 

Director 

US  Army  Materiel  Systems  Analysis  Activity 

Commander 

ATTN:  DRXSY-MP 

1 

US  Army  Dugway  Proving  Ground 

ATTN:  DRXSY-T  (Mr.  Kaste) 

1 

ATTN:  Technical  library,  Docu  Sect 

1 

Aberdeen  Freeing  Ground,  MD  21005 

Dugway,  UT  84022 

Commander 

US  ARMY  TRAINING  4  DOCTRINE  COMMAND 

US  Army  Mfcaflo  Comnmnd 

Redstone  Scientific  Information  Center 

Commandant 

ATTN:  DRSMMtPR  (Documents) 

1 

US  Army  Infantry  School 

Redstone  Arsenal,  AL  35109 

ATTN:  NBC  Division 

1 

Fort  Banning,  GA  31905 

Director  > 

DARCOM  Field  Safety  Activity 

Commandsnt 

ATTN:  DRXOSC 

1 

USAMP4CS/TC4FM 

Charlestown,  IN  47111 

ATTN:  ATZN-CM-CDM 

1 

Fort  ftCcOeOen,  AL  36205 

US  Army  Natick  Reeeerch  4  Development  Commend 

Commander 

ATTN:  DRDNA-VR 

1 

US  Army  Infantry  Center 

ATTN:  DRDNA-VT 

1 

ATTN:  ATSH-CD-MSC 

I 

Natick,  MA  01740 

Fort  Banning,  GA  31905 

US  ARMY  ARMAMENT  RESEARCH  AND 

Commander 

DEVELOPMENT  COMMAND 

US  Army  Infantry  Center 

Directorate  of  Plans  4  Training 

Commander 

ATTN:  ATZB-DPT-PO-NBC 

1 

US  Army  Armament  Reeeerch  and  Development  Conn 

nand 

Fort  Benxdng,  GA  31905 

ATTN:  DRDAR-LCA-L 

1 

ATTN:  DRDAR-LCE 

2 

Commander 

ATTN:  DRDAR-LCEC 

1 

USA  Training  4  Doctrine  Command 

ATTN:  DRDAR-LCU 

1 

ATTN:  ATCD-Z 

1 

ATTN:  DRDAR-LCUCE 

1 

Fort  Monroe,  VA  23651 

ATTN:  DRDAR-PMA  (GJL  Sacco) 

4 

ATTN:  DRDAR-SCA-W 

1 

Commander 

ATTN:  DRDAR-TSS 

2 

USA  Combined  Aims  Center  4  Fort  Leavenworth 

ATTN:  DRCfttCAWfrAM 

1 

ATTN:  ATZLCA-COG 

I 

ATTN:  DRCPMCAWMI 

1 

ATTN:  ATZL-CAM-IM 

1 

Dover,  NJ  07801 

Fort  Leavenworth,  KS  66027 

i  Team 

US  Army  TRADOC  System  Analysis  Activity 

ATTN:  MDAftACW 

1 

ATTN:  ATAA-SL 

1 

AbwdMa  Pnwtag  GraoM,  MD  21010 

Whhi  Sands  MMa  Range,  NM  88002 

la 


DISTRIBUTION  LIST  5  (Contd) 


US  ARMY  TEST  A  EVALUATION  COMMAND 
Commander 

US  Amy  Test  A  Evaluation  Command 
ATTN:  DRSTE-CM-F 
ATTN;  DRSTE-CT-T 
Aberdeen  Proving  Ground,  MD  21005 

DEPARTMENT  OF  THE  NAVY 

Commander 

Naval  Explotive  Ordnance  Dfcpoml  Facility 
ATTN;  Army  Chemical  Officer,  Code  AC-3 
Indian  Head,  MD  20640 

Commander 
Naval  Weapons  Center 

ATTN:  Technical  Library /Code  343 
China  Lake,  CA  93555 

Commanding  Officer 
Naval  Weapons  Support  Center 

ATTN:  Code  5042/Dr.  B.  E.  Douda 
Crane,  IN  47522 


DEPARTMENT  OF  THE  AIR  FORCE 

HQ  Foreipi  Technology  Division  (AFSO 
ATTN:  TQTR 

Wright-Pstterson  AFB,  OH  45433 

HQ  AFLC/LOWMM 
Wrigbt-Patterson  AFB,  OH  45433 

OUTSIDE  AGENCIES 

Battette,  Columbus  Laboratories 
ATTN:  TACTEG 
505  King  Avenue 
Columbus,  OH  4320) 

Toxicology  Information  Center,  WG  1008 
National  Research  Council 
2101  Constitution  Ave,  NW 
Washington,  DC  20418 


US  MARINE  CORPS 

Director,  Development  Center 
Marine  Corps  Development  A  Education  Command 
ATTN:  Fire  Power  DMrion 
Quantico,  VA  22134 


